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Abstract 

A tetra of sets which elements are time series of interbeats has been 
obtained from the databank Physionet-MIT-BIH, corresponding to the fol- 
lowing failures at the humans' heart: Obstructive Sleep Apnea, Congestive 
Heart Failure, and Atrial Fibrillation. Those times series has been analyzed 
statistically using an already known technique based on the Wavelet and 
Hilbert Transforms. That technique has been applied to the time series of 
interbeats for 87 patients, in order to find out the dynamics of the heart. 
The size of the times series varies around 7 to 24 h. while the kind of wavelet 
selected for this study has been any one of: Daubechies, Biortoghonal, and 
Gaussian. The analysis has been done for the complete set of scales ranging 
from: 1-128 heartbeats. Choosing the Biorthogonal wavelet: biorS.l, it is ob- 
served: (a) That the time series hasn't to be cutted in shorter periods, with 
the purpose to obtain the collapsing of the data, (b) An analytical, universal 
behavior of the data, for the first and second diseases, but not for the third. 
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1. Introduction. 

During many years, the study of the Electrocardiogram (ECG) has been 
useful at the time of prognosticating several diseases of the heart. Com- 
monly a Physician uses to observe the ECG and according to its empirical 
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or scientific expertise is able to diagnose an injury at that organ. However, 
that visual study is not always enough to predict or prognostic a sickness 
of the heart because of the high complexity of the EGG. That difficulty 
in the study of that signal comes from the intrinsic aperiodic dynamics of 
the heart and perturbations of that organ due to a set of physiological pa- 
rameters affecting it. Additionally, the noise coming from the electrodes 
due to normal changes of position of a patient could modify dramatically 
the diagnosis of a disease With the purpose of owning an analytical 
method to study the EGG, several statistical techniques have been devel- 
oped many years ago j2|, [sf . One of them is of our interest because unveils 
a hidden behavior of the heart. Such a technique of analysis is useful to 
obtain unusual nonlinear information of the cardiac dynamics of the hu- 
man being js], 0] and another species [sj. The technique follows the steps: 
(i) Select a wavelet (m is the m-th moment, that indicates that 

J^°° x"^ ■ \E'*^'^^(a;) dx = 0) in order to analyze the time series of interbeats: 
Xi = x{ti),i = 1,2, (ii) Select a set of temporal scales: s = l,2,...,So, 
where n G (iii) Apply a Gontinuous Wavelet Transform to the time se- 
ries [el, 0]: Wi = \E'('") ['^J^i'^i ' where Wj are the coefficients 
of the Wavelet Transform, (iv) Apply the Hilbert Transform to the time se- 
ries: hi = Wi+ z Hj, where Hi = H{Wi) 0, Isl, (v) Obtain the amplitudes 
Ai = \/W'l + iff for the elements of the new time series, (vi) Normalize the 
area under the distribution of amplitudes, and (vii) Rescale that distribution 
such as Ai Ar-Pmax and P{Ai) P{Ai)/P^ax, where Pmax is the maximum 
of the normalized distribution. That technique has been successfully applied 
to a set of time series corresponding to healthy persons, resulting after the 
analysis a common Gamma distribution: Gu{x) = b'^^^x'^e'^^ /Tlu+l), where 
b = jy/xo, being Xq the value that maximizes the distribution. The wavelet 
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used is the Gaussian: gn{x) = C'n-r^^ ^ • ^ characteristic parameter for 
the Gamma distribution is u = 1.8 ± 0.1, during the day (6 h.), while during 
the night (6 h.) it is = 1.4 ±0.1. A common Gamma distribution suggests 
that there is an intrinsic dynamics in the heart similar to that of a system 
out of equilibrium jsj. That collective behavior is characteristic also of a 
system that shows a phase transition j^, 0] . It has been found recently, that 
exists a Biorthogonal wavelet so good as the Orthogonal proposed at Ref. [sj, 
that is useful to perform analyzes by the long term (up to the scale of 2^° 
heartbeats) [l3]. If the wavelet biorS.l is chosen a single u = 1.43 ± 0.03 
characterizes the heartbeating of healthy humans during the circadian pe- 
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riod. The robustness of the technique set up by H. E. Stanley et al. together 



the diversity of cardiac data at the pubhc databank Physionet [llj, suggests 
to test both with the goal of learning about the dynamics of the heart in 
patients owners of a failure at that organ. In this work it is studied a triad 



of different sets of patients, suffering of an abnormality at the heart [12 
Obstructive Sleep Apnea (OSA), Congestive Heart Failure (CHF) or Atrial 
Fibrillation (AF), during periods of the order of a day. With the purpose 
to quantify the change in the distribution of normalized amplitudes due to 
the wavelet chosen, that one is selected from any of the three quite different 
families j^: Orthogonal (Daubechies), Biorthogonal, and Non Orthogonal: 
The Gaussian function, which is used by the authors of the wavelet technique 
of statistical analysis. 

2. Obstructive Sleep Apnea. 

This disease consists of the recurrent interruption of the breathing and is 
associated with diverse physiological injuries: drowsiness, high blood pres- 
sure, heart failures, and in general increased mortality. In spite of all these 
symptoms, OSA is quite difficult to detect and usually it is necessary an 
expensive research inside a specialized hospital room. Fortunately there are 
statistical techniques that from the time series of the heartbeat and the 
breath also, are able to find patterns that could avoid those traditionally 



prolonged studies [ij, [15 



3. Congestive Heart Failure. 

This is a dysfunction of the heart that avoids that the flux of blood 



circulates adequately through the organs of the body [12] • Characteristic 
symptoms are: enlarged legs or ankles or difficulty on breathing. Some of 
those symptoms are produced by e.g. narrowed arteries, past heart attack, 
high blood pressure, or heart defects present at birth. There is not reference 
at the hterature of a study of CHF using the statistical technique applied 
here. Interested by answering the same question: Is there a Gamma distri- 
bution that represents the CHF sickness? we did exactly the same procedure 
indicated in previous section but now focused to CHF. 
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4. Atrial Fibrillation. 

This is a disorder of the heart that consists in the quiving of the two 
atria, instead of its normal synchronized beating. The blood isn't pumped 
properly so that the clotes stop its proper fluxing producing a stroke. AF is 
understood also as a stationary electric wave at the ventricles. It is one of 
the first malfunctions of the heart in the world. The intrinsic dynamics of 
the heart in patients suffering of AF is also of our interest. 



5. Analysis of the data. 

The analysis is done for the following databases at Physionet-MIT-BIH 



Databank [llj: Polysomnographic (sip), BIDMC Congestive Heart Failure 
(chfdb). Congestive Heart Failure RR Interval (chf2db), and Atrial Fibrilla- 
tion (afdb). The full set of recordings have been acquired at Boston's Beth 
Israel Hospital (now the Beth Israel Deaconess Medical Center). The data 
is analyzed after choosing the next wavelets Q: Daubechies: dbl, db2, db3, 
Biorthogonal: biorS.l, bior3.3, bior3.5, and Gaussian: gausl, gaus2, gausS, 
following the steps suggested at the Section 1. The range of selected scales 
for the study is 1-128 heartbeats. 

5.1. Obstructive Sleep Apnea. 

A set of 20 OS A patients (sip): men and women between 27 and 60 
years of age, with weights ranging 53-153 Kg. is selected with the purpose 
of studying the statistical behavior of the heartbeat of those sicks [3|. The 
length of the recordings varies in length from 7-10 h., sampled at 100 data 
per second. The RR interbeats where obtained directly from Physionet and 
then applied to them the wavelet analysis. 

5.2. Congestive Heart Failure. 



A set of 15 CHF-1 patients (chfdb, NYHA: III-IV |l3|): 8 females and 
7 males between 22-71 years old is studied. Each one of the recordings is 
around 20 h. in length, sampled at 250 data per second. Another set of 29 



The notation for wavelets is used according to Ref. [6]. For Orthogonal and Non 
Orthogonal wavelets it is used: wavabrN, where small letters indicate the abbreviation 
for the wavelet "wavabr" and N corresponds to the N — 1-th moment. For the Biorthog- 
onal wavelet the notation is: biorNr.Nd, where N^.Nd corresponds to the moments Nr of 
reconstruction and Nd of decomposition. 
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CHF-2 patients (chf2db, NYHA: I-III): 2 females, 8 males, and 21 undefined 
persons between 34-79 years age is analyzed too. That database consists of 
RR time series sampled at 128 Hz. 

5.3. Atrial Fibrillation. 

Here, a set of 23 AF patients has been studied. Individual recordings are 
around 10 h. in duration, sampled at 250 data per second. 

In this study the moments of the wavelets chosen are not considered to 
be higher than 3 because the patients at the hospital are supposed to be 
almost in rest during the day, and that they don't use to follow activities 
that modify their physiological states in a specific way: at least that has 



not been indicated at the databank Physionet [13[. However, during this 
work, in some cases, it is necessary to use moments up to 9 because of the 
complexity of the signals. Analyzes using wavelets in this way have been 



done previously |l7| 



The long term study of the heartbeat of sick people starts applying a filter 
to the time series of RR intervals, according to the algorithm at Ref. 



m 



order to avoid the spurious data. After that step it is followed the procedure 
of wavelet analysis for each of the wavelets and scales already selected. The 
complete analysis of the cardiac time series is done using a software for the 
analysis of ECGs that has been previously described at Ref. [l9j . 



6. Results. 

Results of the analysis for OSA are shown in Figs. 1-4. Good data collapse 
holds for: (a) OSA: Scales: 1-128, Wavelet: bior3.1; (b) CHF-1: Scales: 8- 
128, Wavelet: All of them; (c) CHF-2: Scales: 16-128, Wavelet: Any one 
with moment higher than 3; (d) AF: Scales: 1-128, Wavelet: None. 

The result in (a) shows a quite stable Gamma function through the scales 
of analysis. If it is chosen the wavelet bior3.1 v = 1.23 ± 0.05. This value is 
quite close to the proposed at Ref. in the case of healthy humans monitored 
by night. The intrinsic, universal dynamics of the heartbeat of the patients 
suffering of OSA is now unveiled using the wavelet bior3.1 in the analysis. 
That does not happen selecting any other wavelets already considered in this 
work. 

The result in (b) suggests that the breathing does affect the dynamics of 
the heart at the low scales 1-8, that are in the range of the normal breathing. 
Scales higher than 8 does show data collapse of the times series of CHF-1, for 
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the full set of wavelets applied. As a consequence, due to CHF-1 the intrinsic 
dynamics of the heart is hidden, if it is focused at the scale of breathing. The 
parameter v = 1.56 ±0.07 for representative Gamma distributions, using the 
wavelet biorS.l. On the other hand, result in (c) shows good data collapse 
considering scales of the order of 16 and higher but momenta of the order 
of 3, indicating that the high non stationarity of the EGG corresponding to 
CHF-2 is not easily characterized -perhaps- because a hugger perturbation 
of the breathing on the heartbeat. In this case u = 1.35 ± 0.08. 

The result in (d) indicates that data collapse is not possible applying any 
of the wavelets previously selected. In order to look for the hidden dynamics 
of the heart, higher moments up to 9 in Daubechies and Gaussian wavelets 
are chosen, obtaining a similar answer. Additional Biorthogonal wavelets 
bior2.2, bior3.7, and bior3.9 are unuseful also at the time of looking for 
that universal behavior of the data. That suggests that the dynamics of AF 
patients is so perturbed, that it seems there is not a function representing 
the time series of heartbeatings. 

7. Conclusions. 
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Figure 1: Distribution of the analyzed data corresponding to OSA, after CVAA (shown 
in small black squares here and at the next figures). The parameters of the fit are v = 
1.29 ± 0.04 and = 0.29. 
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Figure 2: Distribution of the analyzed data corresponding to CHF-1, after CVAA. The 
parameters of the fit are v = 1.53 ± 0.06 and = 0.27. 
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Figure 3: Distribution of the analyzed data corresponding to CHF-2, after CVAA. The 
parameters of the fit are v = 1.29 ± 0.06 and = 0.85. 
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Figure 4: Distribution of the analyzed data corresponding to AF, after CVAA. Data 
collapse does not happen. 
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